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\ Abstract. We consider one-dimensional elliptic Ruijsenaars model of type BC\. It is given 

■ by a three-term difference Schrodinger operator L containing 8 coupling constants. We show 

that when all coupling constants are integers, L has meromorphic eigenfunctions expressed 

fT | ■ by a variant of Bethe ansatz. This result generalizes the Bethe ansatz formulas known in 

_ -, ! the Ai-case. 
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1 Introduction 



The quantum Ruijsenaars model p3] in the simplest two-body case reduces to the following 
| difference operator acting on functions of one variable: 

^ ! <r[z) a(z) 

O 

where o~(z) is the Weierstrass a- function, m is the coupling parameter, and T 7 stands for the shift 
operator acting by (T y f)(z) = f(z + 7). This operator, which first appeared in E. Sklyanin's 
work [151116]. can be viewed as a difference version of the Lame operator — d 2 /dz 2 +m(m+l)p(z). 
It was observed by Krichever-Zabrodin [13] and by Felder-Varchenko [U [7] , that in the special 
case of integer coupling parameter the operator ([!]) shares many features with the Lame operator. 
^ . In particular, when m £ Z + they both have meromorphic Bloch eigenfunctions which can be 

given explicitly by a suitable Bethe ansatz. This reflects the well-known fact that the Lame 
operator is finite-gap for integer m (see e.g. [6] for a survey of the finite-gap theory). 

The Lame operator has the following generalization closely related to the Heun's equation: 

3 

H = -d 2 /dz 2 + Y,9 P {g P + l)p(z + u p ), (2) 

p=0 

where u p are the half-periods of p{z). It can be viewed as a £>Ci-generalization of the Lame 
operator. Again, for integer coupling parameters g p this operator is finite-gap, as was discovered 
by Treibich-Verdier Q2], see also [THUS] for the detailed study of ©. 

The operator ([2]) has a multivariable generalization known as the Inozemtsev model [9]. 
A relativistic version of the Inozemtsev model (= BC n version of the Ruijsenaars model) was 
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suggested by J.F. van Diejen [H [5], see also [TUJ [TT]. In the simplest one- variable case it takes 
the form of a three-term difference operator 

L = a(z)T 2 ~< + b(z)T~^ + c{z), (3) 



where 

■J-r a v (z — u v )a v (z + 7 — ui, ) 

-w = II - • w -«(-')■ w 

c(z) is given explicitly in (JSj) below, and the notations are explained at the beginning of the next 
section. 

Therefore, the operator ([3]) should be viewed both as a difference analogue of (|2|) and a BC\- 
version of the Ruijsenaars model ([I]). In the trigonometric limit it coincides with the Askew- 
Wilson difference operator [1]. (Notice that ([3]) contains eight parameters /j, p , fi' p , compared to 
the four in the Askew- Wilson operator and in ([2]).) Therefore, it is natural to expect that (|3]) 
and dT} should have similar properties. This is indeed the case, as we will demonstrate below. 
Our main result says that in the case of integer coupling parameters 

fi p = 2>ym p , n' p = 2jm' p , m p , m' p € Z + (p = 0, . . . , 3), (5) 

the operator ([3]) has meromorphic Bloch eigenfunctions which can be given explicitly via a version 
of Bethe ansatz. Note that our derivation of the Bethe ansatz equations is very elementary: it 
only uses some simple facts about the operator ([3]). 



2 Ruijsenaars operator of type BC\ 
2.1 Preliminaries 

Let cr(z) = a(z; 2lj±, 2t02) denote the Weierstrass a-function with the half-periods uj\, u)2- Recall 
that (j{z) is an entire odd function on the complex plane quasiperiodic with respect to 2oj\, 2co2- 
It will be convenient to introduce the third half period as W3 = —u)\ — LJ2- One has 

a(z + 2uj s ) = -a{z)e 2ris{z+ula \ 8 = 1,2, 3, 

with r] s = ((tv s ), where Q{z) = a'(z)/a(z) denotes the Weierstrass ("-function. Clearly, rji +7/2 + 
773 = 0. There is a relation between r] s and half-periods as follows: 

771W2 - 772^1 = vri/2 . 

It is known that cr(z) has simple zeros at points of the period lattice 21", where T = oj\L + ^2^- 
Likewise, (,(z) has simple poles with residue 1 at those points. Let us introduce the shifted 
versions of a as follows: 

a r (z) = e~ r,rZ a{z + u r )/a(u r ), r = 1, 2, 3. 

These are even functions: a r (—z) = a r (z). The quasiperiodicity properties of a r (z) look as 
follows: 

a r (z + 2lo s ) = (-l) s ^a r {z)e 2ria{z+ul °\ r, s = 1, 2, 3. 
Below we will use the convention that ctq(z) = cr(z) and loq = ijq = 0. Note that 



a(2z) = 2a (z)a 1 (z)a 2 (z)a 3 (z). 
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Let us remark on quasiperiodicity of the coefficients a(z) and b(z) = a(—z) of the operator ([3]). 
Let us use the subscript to indicate the dependence of a(z) = a^{z) and b(z) = bu(z) on the 
parameters \x = {p p ,p' p }. Using the translation properties of a p one checks directly that a(z), 
b(z) given by © have the following covariance with respect to the shift by a half-period u r 
(r = 0,...,3): 

3 3 
Vr E ~Vt E (JH>+Vy) 

a^(z + u) r ) = aji(z)e p =° , b p (z + u r ) = bj L (z)e ?=° , (6) 

with Jlp = fJ-K r (p) and Jl' p = A^ r (p\, where Try is one of the following permutations: 

7r = id, TTi = (01)(23), vr 2 = (02)(13), vr 3 = (03)(12). (7) 

Note that these permutations form an Abelian subgroup of S4, and 7r p (0) = p for p = 0, . . . , 3. 
It follows that Try o ir p = ix q whenever q = Try (p) . 

2.2 Ruijsenaars operator of type BC 1 

Let L be the operator ([5])- (|4"]l with the coefficient c(z) given by 

a p( z ) 

c ( z ) = Z^WpO + 7) - C P (z -7)), C P {z) = —r-; = -r} p + C{z + u p ), (8) 

where c p looks as follows: 
2 3 

Cp = ~^2V) II + ^pW)^Kp(s))- ( 9 ) 
Here the permutations tt„ are the same as in (j7|). 

Remark 1. L is a one-dimensional counterpart of its two- variable version introduced in [2], see 
also [5j 1101 ITT] . Their precise relation is as follows: the -BC2 version in [3] involves one more 
parameter \i attached to the roots e\ =L 62, and it decouples when p = 0. Namely, let Z) be as 
in [4], formulas (3.11), (3.19)-(3.20), (3.23)-(3.24). Then we have, for an appropriate constant 

a, that lim(D — ap -1 ) = L\ + L2, where L±, L2 act in z = Z\ and Z2, respectively, and are given 

jti— »0 

by ©-(H]) and ([H])-© above. Observe that in the special case /Uq = 27m, p' = fx p = p' = 
(p = 1, 2, 3) the operator L reduces to dTJ. 

2.3 Symmetries of L 

It is obvious from the formula dH) for c(z) that c(— z) = c(z). As a result, L is invariant under 
(z <-> — vz). Next, a direct computation shows that the function c(z) = c fJ- (z) is covariant under 
the shifts by half-periods, namely: 

c M (z + uy) = c^(z), where ju p = M*>(p), % = aC( p )- (10) 
Let us write L = L p to indicate dependence on /i. Combining ([6]) and (|10p . we obtain that 

3 

T+ o L M o = e- A ^ o Lji o e A ^, A r = 7 ?r (2 7 )- 1 ^(^ + ^), 

p=0 

where Jl = Tiy(/i) is the same as in ([6]), (|10p . This implies that for any w£T 

roL^o = e ~ A(w)z o L p o e x( - w)z , (11) 

where A (a;) := niAi + 71^2 if a; = n\uj\ + ri2UJ2, and /I in the right-hand side is defined as 
Jl = ir s (p) if uj = to s mod 2r. 
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3 Bet he ansatz 

3 

Providing the coupling constants satisfy ©, put m = ( m p + m 'p) an< ^ consider the following 

p=0 

function ip(z) depending on the parameters t\, . . . , t m , k € C: 

m 

1>(z) =e kz ]Ja(z + t l ). (12) 

i=l 

Let us impose m relations onto these parameters as follows: 

+ 2j 7 ) = ^(lo s - 2j 7 )e 4 ™ (J = 1, • • • , m s ), (13) 

^ s + (2j - 1) 7 ) = - (2j - l) 7 )e^- 2 W (i = 1, . . . , m'J. (14) 

(Here s = 0, ...,3.) We will refer to (|13p - (|14p as the Bethe ansatz equations, or simply the 
Bethe equations. Explicitly, they look as follows: 



in 



n t'T'lfi ^ 4 "* = 1.....™.). (15) 

<7(tj + w s + 2j7) 

y <7(ti + w, + (2j-l) 7 ) w ' ' s; V ; 

Now we can formulate the main result of this paper. 

Theorem 1. Suppose the parameters ti, . . . ,t m ,k satisfy the Bethe equations (I13p ~ (ll4p and the 
conditions ti + tj £ 2u\Z + 2u;2^ for 1 < i ^ j < m. Then the corresponding function tp(z) (fT2|) 
is an eigenfunction of the operator ([3]) -([5]) . 

The proof will be given in the next section. 

Remark 2. To compute the corresponding eigenvalue, one evaluates the expression Lip/ip at 
any suitable point z. For instance, a convenient choice is z = 2 7 mo (provided thq > 0), because 
then the first term in Lip vanishes. 

Remark 3. If some of the coupling parameters m p , m' vanish, then the corresponding sets 
of the Bethe equations are not present in (|13p ~- (ll4p . For example, in the case when the only 
nonzero parameter is too = to, the Bethe equations take the form: 

rp(2jj) = ip(-2jj) (j = 1, . . . , m). 

In that form (seemingly different from |13|, [8] ) they appeared in [21 j . 

3.1 Invariant subspaces 

The idea of the proof of the theorem is that applying L to ip will not destroy the conditions 
(jl3p ~ (|14p . cf. [21 [3]. We begin with two elementary results about a three-term difference operator 
with meromorphic coefficients: 

D = a{z)T 2 ^ + b{z)T- 2 ^ + c(z). 

Suppose that a, b, c are regular at z € 2 7 Z, apart from z = where a, b have simple poles. 
Furthermore, suppose that 

res 2=0 (a + 6) = (17) 
and that for some to £ Z + the following is true: 

a(2 7 TO) = 0, a(2j 7 ) = 6(-2n), c{2j 1 ) = c(-2j 7 ) (j = ±1, . . . , ±m). (18) 
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Lemma 1 (cf. [2], Lemma 2.2). Let D be as above and define Q m as the space of meromorphic 
functions f(z) which are regular at all points z 6 27Z and satisfy the conditions f(2jj) = 
/(-2j7) for all j = 1, . . . , m. Then D(Q m ) C Q m . 

Proof. For D' = aT 2 ^ + bT~ 2 "' this is precisely Lemma 2.2 from [2], thus D'(Q m ) C Q m . On 
the other hand, the conditions on c in (|18|) imply trivially that cQ m C Q m . ■ 

Corollary 1. Suppose that instead of (|17p . (|18p u>e &nou> £/iai D is invariant under (z <-> — z) 
and i/iai a(2jm) = 0. TTien D{Q rn ) C Q m . 

Proof. Indeed, in that case we know that = a(— z) and c(z) = c(—z) identically. This 
implies the conditions (|17 p - (|18p . ■ 

For the next lemma, we assume that: (1) a is regular at z € 7 + 27Z apart from a simple 
pole at z = —7; (2) b is regular at z 6 7 + 27Z apart from a simple pole at z = 7; (3) c is 
regular at z € 7 + 27Z apart from simple poles at z = ±7. Also, suppose that 

res 2= _ 7 (a + c) = res z=1 (b + c) = 0, res 2= _ 7 a = res z=7 6, (19) 
(a + 6 + c)| 2= _ 7 = (a + b + c)\ z=1 . (20) 

(The last condition makes sense because f)19|) implies that a + 6 + c is regular at z = ±7.) In 
addition to that , assume that for some m € Z+ the following is true: 

o((2m - 1)7) = 0, o((2j - 1)7) = 6((-2j + 1) 7 ) for j = ±1, . . . , ±m, (21) 
c((2j + 1)7) = c(-(2i + 1)7) fori = l,...,m-l. (22) 

Lemma 2 (cf. [2], Lemma 2.3). For D as above, define Q' m as the space of meromorphic 
functions f{z) which are regular at z E 7 + 27Z and satisfy the conditions f({2j — 1)7) = 
f((-2j + 1)7) for j = 1, . . . , m. Then D{Q' m ) C Q' m . 

Proof. For D' = a(T 2 ^ - 1) + 6(T- 2 ^ - 1) the proof of the inclusion D'(Q' m ) C Q' m follows the 
proof of Lemma 2.3 in [2] . On the other hand, the conditions on a, b, c imply that d := c + a + b 
belongs to Q' m . Thus, c'Q' m C Q' m . Therefore, the operator D = D' +{a+b+c) preserves Q' m . ■ 



Corollary 2. The lemma above remains valid after replacing (|20p - (|22p by the invariance of D 
under (z *-> —z) and the condition that a((2m — 1)7) = 0. 

Proof. Indeed, the conditions (120[) — (j22f) follow easily from the fact that b{z) = a(—z) and 
c(z) = c(—z). ■ 

Let us apply these facts to the Ruijsenaars operator (|3|) with integer coupling parameters ([6]). 
Below we always assume that the step 7 is irrational, i.e. 7 ^ Q <g>z T = Qloi + Qu>2- We proceed 
by defining Q as the space of entire functions ip( z ) satisfying the following conditions for every 

weu s + 2r (s = 0, ... ,3): 

1>(u + 2j 7 ) = - 2n)e^ mr >^ (j = 1, . . . , m s ), (23) 
^ + (2j - 1)7) = ^{u - (2j - l)7) e W-2)7m,(<,) (j = 1, . . . , m ' s ). (24) 

3 

Here m stands as before for m = ^ (m p + Tii), and the constant is defined for lv = 

p=0 

niui + n 2 uJ2 as r)(u) = n\r\\ + n 2 ri 2 . 

Proposition 1. For integer coupling parameters ([6]) the Ruijsenaars operator ([3]) preserves the 
space Q of entire functions with the properties ([23]) - ([24"]) : L(Q) C Q. 
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Proof. First, by applying Corollaries [H [2] to the Ruijsenaars operator, we obtain that L pre- 
serves the spaces Q mo and Q' , (in the notations of Lemmas ([T]), Note that in doing so, 
we only have to check the vanishing of a{z) as required in Corollaries [H [21 and the conditions 
on the residues (|19p . This is where the formula ([9|) becomes crucial. Finally, in order to show 
that L preserves similar conditions at other points w£f, one applies the formula (fTT|) . ■ 



Next, given ol\,oli £ C and m £ Z, let us write J-m ,a2 for the space of meromorphic func- 
tions ip(z) having the following quasiperiodicity properties: 

ib{z + 2u s ) = e m ^ z+a °ip(z), a = 1,2. 

Entire functions in J-^ ,a2 (m > 0) are known as theta- functions of order m (with characteristics), 

each of them being a constant multiple of (|12p . for appropriate t\, . . . , t m , k. 

Now, a simple check shows that in the case ([6]) the Ruijsenaars operator ([3]) preserves these 

3 

spaces corresponding tom = ^ (m p + m' p ): 

L{T^)^T^ a \ Vai,a 2 . 

Combining this with Proposition [TJ we conclude that L preserves the space of theta-functions 
of order m satisfying the conditions (|23p ~ ([24 



L(C- a2 nQ)Cf^nQ, Vai,a 2 . (25) 

Proof of the Theorem [TJ Take a solution (t\, . . . , t m , k) to the Bethe equations and the cor- 
responding function ip (]12p . Clearly, ip belongs to Tm ,a2 for some ol\,oli. The Bethe equations 
give the conditions (|23p - (|24p only for uj = u s , but the rest follows from the translation properties 
of ip. Thus, ip belongs to the space J^ 1 '" 2 fl Q. By (|25]) . ip := Lib also belongs to this space. 
Now we use the following fact (whose proof will be given below): 

Lemma 3. For any two functions ip,ip & Fm' a2 H Q, their quotient ip/ip is an even elliptic 
function, i.e. it belongs to C(p(z)). 

By this lemma, if ib/ip is not a constant, then its poles must be invariant under z *— > —z, thus 
there exist at least two of t\, . . . , t m such that their sum belongs to 2r. ■ 

Proof of the lemma. Take any two functions ip, ip in J-m'" 2 HQ and put / := ip/ip. Note that 
/ is an elliptic function of degree < m (because its denominator and numerator have m zeros 
in the fundamental region). Let us label m pairs of points uj s ± 2jj, lu s ± (2j — 1)7 as P, with 
I = 1, . . . , m, then the properties of ip, ip imply that / satisfies the conditions 

f(P+) = f(Pf), l = l,...,m. (26) 

We may assume that / is regular in at least one of the half-periods uj s , otherwise switch to 
1/f = ip/ip. Let us anti-symmetrize / to get g{z) := f(z) — f(—z), which will be odd elliptic, of 
degree < 2m. It is clear that g also satisfies the conditions (|26p . At the same time, it is anti- 
symmetric under any of the transformations z 1— > 2u s — z [s = 0, . . . , 3). Altogether this implies 
that g must vanish at each of the 2m points P^. Finally, it must vanish at one of the half-periods 
(where / was regular). So g has > 2m + 1 > deg(g) zeros, hence g = 0, f(z) = f(—z), and we 
are done. 

The above argument, however, would not work if one or both of the functions ip, ip vanish at 
some of the points PrK Indeed, then we cannot claim that / is regular at those points, so some 
of the conditions (f26|) would not hold for /. In that case, we can argue as follows. Let ip\ denote 
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the linear combination ip\ = ip + Xtp. Then ip\, tp^ for generic A, fi will have zero of the same 
multiplicity at any given point P^. Thus, choosing A, fi appropriately, we can always achieve 
that ipx/ip/j, 7^ 0,oo at every of these 2m points. Let r be the number of those pairs (P t , Pf) 
where ip\, ip^ vanish. Then we have that their ratio / := VOv/^/V still satisfies the conditions (|26p 
at the remaining m — r pairs of points and has degree < m — r due to the cancelation of the 
zeros in the denominator and numerator of /. Thus, the previous argument applies and gives 
that / is even. Therefore, tp/tp is even. ■ 

3.2 Continuous limit 

As remarked in [3], the operator ([3]) with the coupling parameters ([6]) in the continuous limit 
7 — > turns into the £?Ci-version (|2j) of the Lame operator: 

3 

L = const + j 2 w(z) o H o w~ l (z) + o(7 2 ), where w{z) = (o p {z)) 9p , (27) 

p=0 

and the coupling parameters g p are given by g p := m p + m' . 

3 

To formulate a Bethe ansatz for the operator ([2]), we put m = Yl g p and let ip{z) = 

ift(z; k, ti, . . . , t m ) be the function (|12p . Let us impose the following m relations on the pa- 
rameters k,t±, . . . , t m : 



d 2 i- 



1 



dz 2 i 



T (i>(z)e- m ^) 



for j = 1,... ,g s and s = 0, . . . , 3. (28) 



Theorem 2. Suppose the parameters t±, . . . , t m , k satisfy the Bethe equations (|28|) and the con- 
ditions ti + tj £ 2oj\L + 2^2^ for 1 < i ^ j < m. Then the function w~ l {z)ip{z) given by (|12p 
with w as in (|27p , is an eigenfunction of the operator (|2|) . 

This theorem is proved analogously to Theorem [TJ 

Example 1. Let g p = for p = 1,2,3 and go = m. Then the operator ([2]) becomes the 
Lame operator —d 2 /dz 2 + m(m + l)p(z). Its eigenfunctions have the form ip{z)a~ m {z) with 

m 

ip(z) = e kz Y[ °~{ z + tj)- The Bethe ansatz equations (|28"|) for k,t\, . . . , t m in this case reduce to: 
i=i 

^_^(0) = for j = l,...,m. (29) 

We should note that this particular form of Bethe equations differs from the classical result by 
Hermite [20| . For instance, in Hermite's equations one discards the points with ti = tj mod 2r, 
while in Theorem [2] we discard the points with ti = —tj mod 2r. Thus, comparing Theorem [2] 
with the Hermite's result, we conclude that (f29j) must be equivalent to Hermite's equations [20] 
provided ti ± tj ^ 2r for i 7^ j. The same remark applies to the equations (|28p when compared 
to the Bethe ansatz in, e.g., [T8]. 



3.3 Spectral curve 

Let us say few words about the structure of the solution set X C C m x C to the Bethe equations 
([T3|) - (fTT|l . We will skip the details, since the considerations here are parallel to those in [HI El 

El EE]. 
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First, using the properties of cr(z), one observes that ip(z) acquires a constant factor under 
the transformations 

(ii, ■ ■ ■ , t m , k) i— > (ti, . . . ,tj + 2w s , ... ,t m ,k - 2i] s ) (s = 1,2). 

As a result, X is invariant under these transformations. Also, multiplying ip(z) by e mz h does 
not affect the Bethe equations, because such an exponential factor is (anti)periodic under the 
shifting of z by multiples of 7. Therefore, A is invariant under the shifts of k by iri/'y: 

(ii, . . . , t m , fc) 1 ^ (ti, . . . , t m , fc + 71-1/7). 

Finally, ^ does not change under permutation of t± , . . . , t m , so A is invariant under such per- 
mutations. 

Let X denote the quotient of X by the group generated by all of the above transformations. 
Explicitly, let b s j(ti, . . . , t m ) and b' s j(ti, . . . , t m ) denote the left-hand side of equations (fT5l) 
and (fT6|) . (Here s = 0, . . . , 3 and j = 1, . . . , m s or j = 1, . . . , m' s , respectively.) Introduce the 
variable q := e 2lk . Then X is described by the equations 

b sd (t 1 ,...,t m )=q 2j , b' sJ (t 1 ,...,t m )=q 2 ^ 1 . (30) 

Excluding the q- variable from the equations (|30p , we may think of X as an algebraic subvariety 
in the symmetric product S m £ of m copies of the elliptic curve £ = C/2T where T = 7Luj\ + Zuj2- 
(See, however, Remark 0] below.) Counting the number of equations, we conclude that every 
irreducible component of A has dimension > 1. Since we are interested (cf. Theorem [T]) in those 
points (ti, . . . , t m ) of A where U + tj ^ 2r, we should restrict ourselves to the open part Y C A, 
lying in 

S m 8 \ U 4<i {tj + tj = mod2r}. 

We need to show that Y is nonempty and one-dimensional. To this end, one easily observes from 
the equations (fT5|) - (fT6]) that the closure Y of Y in S m £ contains the points P + = {P^~ , ■ ■ ■ , P£) 
and P~ = (Pi , ■ ■ ■ ,Pm)i m the notations of the proof of the Lemma El These 'infinite' points 
correspond to q — ► 0, 00 in (|30p . Similarly to [7], lemma 3.2, one shows that near P ± the 
variety Y looks like a smooth curve, with q ±l being a local parameter. One can show that Y 
is an irreducible, projective curve, and it should be regarded as the 'spectral curve' for the 
operator L. For every (t\,...,t m ) S Y \ {P + ,P~}, the corresponding value of q = e 27fc is 
determined from (|3U|) . and the corresponding ip(z) is unique, up to a factor of the form e mNz h , 
We have Lip = eip, with the eigenvalue e being a single-valued function on Y which has two 
simple poles at P^. There is an involution v on Y, which sends (t±, . . . , t m ) to (— ti, . . . , — t m ) 
and the corresponding tfj(z) to ip(— z); note that v(P + ) = P~. The function e is z^-invariant, 
and takes each its value exactly twice on Y. It is straightforward to compute the asymptotics 
of e and ip near P ± . Finally, for generic value of e, the eigenspace of meromorphic functions 
{/ : Lf = ef} is spanned by the corresponding ip(z), ip(—z) over the field K of 27-periodic 
meromorphic functions of z. 

Remark 4. Note that in the case when all m' s = 0, the second set (fT4"|) of the Bethe equations 
is absent, thus a shift k 1— ► k + |^ is also allowed. In that case the coefficient c(z) (JSj) vanishes, 

7TZZ 

so L has two terms only, and it is easy to see that the transformation ijj <— > e 2 i ip changes the 
sign of the eigenvalue e. As a result, the subvariety of S m £ which was obtained by excluding q 
from ([30]), will be a quotient of A by Z 2 , rather than A itself (cf. H3 EI] ) . 
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